Spectrum of the Andreev Billiard and Giant Fluctuations of the Ehrenfest Time 
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The density of states in the semiclassical Andreev billiard is theoretically studied and shown 
to be determined by the fluctuations of the classical Lyapunov exponent A. The rare trajectories 
with a small value of A give rise to an anomalous increase of the Ehrenfest time te ~ | ln?i|/A and, 
consequently, to the appearance of Andreev levels with small excitation energy. The gap in spectrum 
is obtained and fluctuations of the value of the gap due to different positions of superconducting 
lead are considered. 
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Introduction. — The density of states in a metallic is- 
land coupled to a superconductor is modified due to the 
proximity effect [Ij . The changes are most pronounced in 
the vicinity of Fermi energy, where there opens a gap in 
the spectrum of excitations. A ballistic chaotic normal 
region coupled to a superconductor via the small con- 
striction (NS interface) is called the Andreev billiard . 
The spectrum of such billiards was calculated a decade 
ago P| assuming the random matrix description of the 
quantum dynamics. This approximation is valid if the 
number of channels N supported by the NS interface is 
small. In spite of large efforts 0,1111111111113113, 
there is currently no reliable calculation of the spectrum 
in the most interesting case of the semiclassical Andreev 
billiard, when the number of open superconducting chan- 
nels scales with the Planck constant as N ~ 1/H. 

The properties of Andreev billiards are governed by 
the Andreev reflection, when an electron trajectory is re- 
traced by the hole that is produced upon absorption of a 
Cooper pair at the NS interface. At the Fermi energy Ep, 
the classical dynamics of the hole is the time reverse of 
the electron dynamics, so that the motion is strictly pe- 
riodic. This periodicity, however, comes in conflict with 
a quantum-mechanical evolution, since each Andreev re- 
flection adds an extra phase ir/2 to the electron- hole wave 
function. This phase is compensated by the difference of 
classical actions of electron (Ep + e) and hole (Ep — e) 
along the trajectory, where e is the excitation energy. 
The longer the interval t is between Andreev reflections, 
the smaller energy suffices to produce the missing phase 
e = h%/2t. On the other hand, the probability of particle 
trajectory to not touch the NS-interface for a long time 
become exponentially small ~ e~*'* D , with tp, being a 
dwell time. This leads to the prediction of an exponen- 
tial suppression of the density of states H 



p(e) ps N(Tnr/e 2 t D ) exp (-hn/2et D ) 



(1) 



However, the density of states at small excitation ener- 
gies, corresponding to large times, is not captured by this 
formula. A new time scale responsible for the spectrum 
at low energies £| is the Ehrenfest time Tp ~ |ln7i|/A, 
where A is a Lyapunov exponent in the normal bil- 
liard. For longer times the initial quantum wave packet 



AxAp ~ h acquires a macroscopic size due to exponen- 
tial ~ e xt divergency of trajectories. This invalidates the 
trajectory based derivation that led to Q. It is believed 
that, below certain energy £ gap ~ h/rp, the density of 
Andreev states vanishes exactly, p(e < £ gap ) = 0, but the 
magnitude of the gap and the mechanism of its formation 
remained a subject of controversial discussion [rjl^Hoj. 

At finite times, the value of the Lyapunov exponent A 
depends on the specific trajectory [l|, leading to fluc- 
tuations of the Ehrenfest time. Both the Andreev spec- 
trum [isiiHsiinmiii, and the quantum to cias_ 

sical crossover in ballistic transport El El El El El 
attracted a great deal of interest recently, but the role 
of fluctuations of the Ehrenfest time was never investi- 
gated. In this Letter we show how the low energy den- 
sity of states is determined by the large Ehrenfest time 
fluctuations and solve the long standing problem of the 
Andreev gap. 

The distribution of finite time Lyapunov exponents is 
parameterized as (here To is, e.g., an averaged time be- 
tween bounces at the walls of normal billiard) 



P(X ) t)=T exp[tF t (X)}. 



(2) 



In the case of chaos the limit i*t>T (A) = F(X) ex- 
ists, with F(X) specific for the dynamical model [2fj| . 
The function F(X) has a maximum at A = Ao, which 
is the conventional self-averaging Lyapunov exponent 
[-F(Ao) = 0]. Since all small values of A > are present 
in the distribution Pjl. one may always find (rare) tra- 
jectories with any large value of the Ehrenfest time. To 
build a semiclassical eigcnfunction of 2-dimensional bil- 
liard, however, one needs a family of trajectories, all hav- 
ing the same interval t between the Andreev reflections. 
The explicit construction of discrete Andreev levels from 
the tube of trajectories (whose transverse Poincare sec- 
tion is quantized via the Bohr-Sommerfeld rule) is pre- 
sented in Ref. . The gap in the spectrum is determined 
by the largest time for which the number of trajectories 
is enough to form at least one eigenstate. Quantitative 
counting of the number of trajectories is done with the 
use of the concept of transmission band [TH IT^ | . This 
lead us to the expressions for the gap and the density of 
states [(fT5)l and ljH} below] depending on F(X), tp,, and 
a number N of open channels in the NS interface. 
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FIG. 1: Distribution of Andreev levels with e = 7i7r/12ro and 
e = hn /14ro over transmission bands of different size, corre- 
sponding, respectively, to the 6th and 7th iterations of the 
standard map averaged over 200 positions of the NS inteface. 
Dwell time to = 10; kicking strength K — 10. The dashed 
lines show the theoretical prediction 110L with the function 
F found numerically for 5 and 6 iterations. For M = 10 r , the 
value A = would correspond to Nb ~ 10 5 in Eq. @. 

Stroboscopic model. — All of the essential features of a 
generic Andreev billiard are captured by the stroboscopic 
model of the Andreev billiard, which was developed in 
Ref. in order to use the advantages of an open kicked 
rotator |2l| for large scale quantum-mechanical simula- 
tions. Here we present a semiclassical solution of this 
model. First, the quantum kicked rotator, a counterpart 
of the classical standard map , 

p n+ i =p n +KI sin 6>„/t , 6 n+ i = 9 n +T p n+ i/I , (3) 
is defined by the Floquet operator 



,. , Mo d 2 \ ( ,KI cos( 



(4) 



Here Iq,tq, and K are the moment of inertia, the inter- 
val between kicks, and the kicking strength, respectively. 
Next, introduce the dimcnsionlcss Planck constant h e g = 
Htq/Iq. If h e ff = 2ir/M with integer M, the coordinate 
and the momentum p = —ih c gd/d9 take discrete values 
9 k = 2irk/M, p m = 2Trm / A'l , k,m = 1,2,. ..M. The 
Floquet operator now becomes M x M matrix. 

The electron and the hole components of the wave 
function of Andreev kicked rotator span over the doubled 
2M-sites Hilbert space with their evolution given by the 
normal (U) or conjugated (U*) Floquet operators. The 
electron is converted into the hole by reflection at the N- 
channel superconducting lead, attached at 9\ < 9 < 92 
(02 — ®\ = UN). This is done with the help of projection 
matrix Q those N only nonzero elements are Qk,k = lj 
for 9\ < hk < 02. Andreev levels are found from 



- e leTo/h ib U-( ^ ~ ® U ~ 
~ 6 _ iQU (1 



iQU* 
-Q)U* 



(5) 



The classical limit corresponds to M, N — > oo, while the 
dwell time to = tqM/N is fixed. A classical particle 
at any time t — utq has a definite position either inside 
the normal region or at the interface. A semiclassical 
quantization of the map JSJ requires a construction of 
the quantum states having a similar property. A formal 



description of the wave packet <fi, which is injected from 
the superconductor, stays inside the billiard for n — 1 
kicks, and then hits the NS interface is given by (0 < 

771 < 77.) 



Q<j> = cj) , QU m cj) = , QU n cf) = U n <j>. 



(6) 



Provided such a solution is found, one easily builds 2n 
eigenfunctions with the eigenvalues (r = 0, n — 1) 



±hn- 



2r + l 

277To 
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£ n o with the largest possible 77 constitutes a gap. Below 
we always consider only the levels with r = 0. Equa- 
tions © and J7J) are the analog of adiabatic quantization, 
developed in ref. [|| for the generic Andreev billiard. 

Strictly speaking, Eqs. JSJ for any 77 have no solutions. 
However, for h — > where exists « N linearly indepen- 
dent wave packets satisfying © with practically any de- 
sired accuracy, <j^QU m 4> ~ e~ x / h . Finding the number 
of such solutions reduces to the calculation of certain 
phase-space areas called transmission bands flU - 

Transmission bands. — We call the transmission band 
a simply connected part of the phase-space area of the 
NS interface, 9\ < 9 < 9 2 ,0 < p c ff < 2tt, each point 9,p 
of which visits the interface at the n-th iteration of the 
map (and do not visit it earlier). The image of the stripe 
9\ < 9 < 92 after 77 iterations is another long and narrow 
(curv ed) strip e of a width ~ e~ xt N/M (see examples in 
Refs. [HIII3)- Phase-space overlaps of the NS interface 
with its image, the transmission bands, are the areas with 
approximately the shape of parallelogram whose long and 
short sides have a length - N/M and ~ e~ xt N/M. The 
number Nb of (families of) Andreev levels J7J supported 
by a single transmission band is calculated as its area 
divided by 27r/7 ff, 



N b « (N 2 /M)e 
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(8) 



The total number of levels composed from trajectories 
having the time t = utq between hitting the NS interface, 
whose Lyapunov exponent falls in the interval dX, is @ 



dN t = NP(X,t)e- t/tD (T /t D )d\. 



(9) 



These levels originate from the many transmission bands 
of various size: N t — ^ iVj. We may use (JSJl to express 
A through Nb and to find the distribution of levels over 
the sizes of the bands 



dNt N 2 T /.p/l, N 2 

= exp < tbt — m 

dliiN b Mt \ \t MN b 



The distribution dNt/dhiNb, found numerically for the 
model is shown by the histogram in Fig. I for times 
t = 610 and t — 7tq. The choice of 7j e ff = 27r x I0~ 7 
introduces a quantum-mechanical scale in classical area 
counting. Direct quantum-mechanical calculation of en- 
ergy levels in such an Andreev billiard would require di- 
agonalization of the matrix of the size M = I0 7 , which 
is beyond the reach for existing computers. 
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FIG. 2: The number of Andreev levels vs logarithm of the 
area of transmission band dNt/dlnN), for the stroboscopic 
model with K — 11, h c g = 2n x 1CP 7 . The superconduct- 
ing lead is attached at tv/10 < 8 < 37r/10. Histograms for 
different number of kicks n = 2 — 16 are offset vertically and 
multiplied by 2™. For n > 6, theoretical distribution HOH is 
also shown (dashed line). The gap amounts to e™, = 7i7r/3uTo 
instead of e gap = %-k/IAto expected from Refs. (b. ItI Ho|l . 

The dashed lines in Fig 1 show the predicted distribu- 
tion [Eq. (|10|) ] with the numerically found function -Ft (A) 
for times t — 5tq,6tq. The number of iterations in the 
theoretical formula is reduced by 1 because the effective 
strong (~ K) stretching of the phase-space image of the 
NS interface starts only from the second iteration. 

Even though the trajectories with A ~ are present 
for both times shown in Fig 1, the largest transmission 
bands expected from eq. © at the smallest values of the 
Lyapunov exponent are absent [2j|. Also, the area of the 
largest existing transmission band is smaller for a larger 
time. This is the actual reason for the formation of the 
gap in the excitation spectrum. 

The gap in spectrum is given (|7J) by the longest time 
for which there exist the transmission bands with Nt > 1 



£ga P = hn/2t TI 



(11) 



To illustrate the mechanism of the creation of the gap 
we show in Fig. 2 the distribution dNt/dlnNi, for times 
</to = 2— 16 (shifted vertically). Counting the small ar- 
eas numerically is a very time consuming-procedure |25| . 
Starting from the 10th iteration more than half of An- 
dreev levels are missing due to the errors in area count- 
ing. Still, the results for the largest transmission bands 
are trustable for all times presented on the figure. The 
existing estimates of the Andreev gap 0, 0, El neglected 
the fluctuations of the Lyapunov exponents (shown by 
narrow peaks connected by a line). The right solid line 
in Fig. 2 shows the (asymptotic) time dependence of the 
largest transmission band area. The ratio of the slopes 
of two lines give the suppression of the Andreev gap due 
to the Ehrcnfcst time fluctuations. 

Density of Andreev states. — The area of transmission 
band equation (JSJ) increases exponentially for small val- 
ues of A, while the number of trajectories with small A is 



exponentially small ©. These two conflicting effects al- 
low one to find the value of the Lyapunov exponent lead- 
ing to the largest, for a given time, transmission band. 
Equating via (|1U|I the total number of expected levels 
with A < A c to the area of a single band I© with A = A c 
we find 



F(X C ) = -X c + t 



D 



The value of the gap now is 



(12) 



(13) 



The value of A c depends on details of the specific model. 
It was shown in Ref. [24| that the derivative ^'(A) has a 
maximum at A = and that F'(0) < 1. Since the only 
maximum of the function F is F(Xq) = we obtain 



A c < A /2. 



(14) 



;ev ga p is ; 

twice smaller than predicted previously 0, HI lldf . 

Semiclassical methods may be applied for a description 
of eigenstates constructed from the trajectories with the 
Lyapunov exponents 



A < X m ^ = t- 1 \n(Nr /t D ). 



(15) 



For A > A max , the number of levels per transmission band 
became less than 1. The number of Andreev levels asso- 
ciated with the time t may now be found by integration 
of Eq. (O over the interval A c < A < A max . If A max > Ao 
this allows one to recover the known result eq. Q , which 
is now valid for eln(NTo/trj) > hirXo. The smaller ener- 
gies may exist only due to A < X ma x < Ao, for which we 
found a novel form of the density 



P(e) 



Nrn 



et D F'{n) 



exp 



hit , „ hir 
2S F{ ^ ~ 2^ 



(16) 



where k = (e/Tiir) ln(A^To/t_o). Equation (|16|l is valid for 
A c < eln(NTo/t D )/hir < A . 

Formulas (|13fl and (|16f) are the main results of this 
Letter. They describe not only the stroboscopic model, 
but any Andreev billiard with chaotic dynamics coupled 
to a superconductor through the A^-channel lead. (In 
this case, To may be replaced by the averaged time be- 
tween bouncing of the billiard walls. The precise value 
of To is not important since ~ h^ 1 3> tp/To-). For 
the model 10 the semiclassical density consists of a se- 
ries of (5-function peaks at e = e nr Q, and Eq. i|16|) de- 
scribes the smoothed envelope of this distribution. Such 
peaks in the Andreev spectrum were seen in simulations 
of Refs. HEll. 

In Eq. i|T3|) , we found the averaged value of the super- 
conducting gap. Fluctuations of £ gap are caused by the 
variations in the position of NS interface. These fluctu- 
ations result from the fluctuations of the area of largest 
transmission band at a given time, which is shown in 
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FIG. 3: The area of the largest transmission band (in units 
of h — 2ty x 10 -7 ) for 9,10,11, and 12 iterations as a function 
of the position of superconducting contact. Because of the 
symmetries 9 — » — 9,6 — n ^ n — 9, only a half of the [0, 2ir] 
range of variation of the argument is shown. 

Fig. 3. Since the area of (even the largest) transmission 
band decreases exponentially with time while the num- 
ber of the bands exponentially increases, we expected 
that the correlation length for fluctuation of the largest 
area should also decrease with time. No such increased 
sensitivity is seen in Fig. 3. 

Discussion. — Among the total number of M Andreev 
levels, M q levels could not be described semiclassically, 

M q « M 1 '^ (1 + \nN/Xt D ) fj^ 5 > 1. (17) 
These levels originate from the phase-space area of the 



NS interface covered by small transmission bands, N\, < 
1. References 0,EJ predict the vanishing of the density 
of these levels for e < hwXo/2\nN. Although these pa- 
pers do not provide a rigorous calculation of the Andreev 
spectrum, we refer to these results as an indication that 
the levels missing in the adiabatic quantization Q do not 
change the low energy density of states Ijl6(l found in this 
Letter. 

The density of Andreev levels (jlfij) and the Andreev 
gap (|13(l . which we found in this Letter, depend on 
the entire distribution of finite time Lyapunov expo- 
nents P{X,t) [or F(X) J2J], not on the most prob- 
able Lyapunov exponent Ao, as was expected previ- 

ousiy a a saB ana. 

The variations of the Lyapunov exponent leading to 
the fluctuations of Ehrenfest time considered in this Let- 
ter are of the order ~ Ao. However, the corresponding 
variations of the extent of the divergency of trajectories 
become exponentially enhanced ~ e At . So these are in- 
deed the giant fluctuations. The presented results are 
supported by the classical numerical simulations. Ver- 
ification of our findings in a real experiment or in a 
quantum-mechanical simulations 26] remains a challeng- 
ing problem. 
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